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Abstract 


The mechanism of methane dissociation on an Rh-decorated Ni(100) surface has been investigated Using density 

functional theory. The study includes the determination of the most stable adsorbate/adsorbent configurations of the 

species associated with subsequent reactions and generating the energy surface for CH, dissociation process. The Rh- 

decorated Ni(100) surface was found to be more favorable for the process than the NiRh(111) configuration, mainly due 

to lower the activation energy of CH decomposition reaction by 48.5%, leading to a higher conversion of CH, to carbon 

and hydrogen 
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Introduction 


The production of hydrogen and carbon via the catalytic cracking of methane has attracted the attention of many researchers (Abbas 
and Daud, 2010; Abbas and Wan Daud, 2010; Amin et al., 2012; Amin et al., 2011; Ashik et al., 2017; Fan et al., 2012; He et al., 2018; 
Italiano et al., 2010; Li et al., 2018; Serrano et al., 2010; Tezel et al., 2019; Zhang et al., 2017; Zou et al., 2019). Methane is hydrogen- 
rich with a lower carbon supply rate, which can reduce the process of graphene growth on the catalyst (Yang and Chen, 1989; Zhu et 
al., 2007). Transition metals, notably, Ni-based catalysts have shown high efficiency for methane dehydrogenation reaction (Amin et 
al., 2012; Fan et al., 2012; He et al., 2015; Li et al., 2015b; Liu et al., 2011; Shadravan et al., 2018; Wang et al., 2006; Zhu et al., 2007). 
The mechanism of CH, dissociation has been investigated both experimentally and theoretically. Using high-resolution electron energy 
loss spectroscopy (HREELS) (Lee et al., 1987; Lee et al., 1985) and X-ray photoelectron spectroscopy (XPS) techniques (Kaminsky et 
al., 1986; Yang et al., 1995), chemisorbed CH3, CH}, and CH species have been observed on the Ni surface indicating that CH, is 
transforming into C and H, through sequential CH(x:4-1) dehydrogenation reactions. The mechanism of methane cracking appears to 


have the following steps: 


CHa g) + 2 *= CH3(qas) + Fads) (1) 
CH3adsy +*= CHzçaas) + Hoaas) (2) 
CHoadsy +*= CH (ads) + Hoaas) (3) 
CH ads) +*= Caas) + A(aas) (4) 


where the asterisk (*) denotes an active site on the surface, and (ads) means adsorbed species on the surface (Li et al., 2015b). 
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Bimetallic Ni-based catalysts were recently used to improve catalyst stability and activity. Various studies have shown that bimetallic 
Ni-based catalysts are more durable and lead to higher CH, conversion levels (Fan et al., 2012; He et al., 2015; He et al., 2018; Horvath 
et al., 2011; Jia et al., 2020; Li et al., 2013; Liu et al., 2011; Tezel et al., 2019; Theofanidis et al., 2020; Tian et al., 2018). Rhodium is 
known as a suitable catalyst for the hydrogen evolution reaction (HER). However, due to its high cost, many studies have been carried 
out to develop a catalyst characterized by high performance with low Rh content (Huang et al., 2017; Nguyen and Choi, 2019; Nguyen 
et al., 2019). Li et al (Li et al., 2013) have studied the dissociation of CH, on the NiRh(111) surface and compared it with various 
bimetallic Ni-based catalysts using DFT calculations. They have reported that NiRh have higher catalytic activity for CH, dissociation 
than other catalysts. The effect of the orientation of the catalyst surface on methane decomposition was investigated. Yang and Chen 
have shown that the Ni(100) and Ni(110) are among the gas/metal interface, whereas Ni(111) and Ni(311) belong to the graphite/metal 
interface. In their study, they concluded that the activation energy required for CH, dissociation on Ni(100) was lower than that on 
Ni(111) (Zhu et al., 2007). In this study the (100) orientation is incorporated with doping Rh with Ni to investigate the mechanism of 
CH, cracking process. The aim is to show that the Rh-Ni(100) surface is more adroit for the abovementioned process compared to the 
NiRh(111) configuration, particularly, for the CH dissociation step. 


1 Materials and Methods 
1.1 Computational Method 


Density functional theory calculations were carried out to investigate the mechanism of CH, dissociation on the Rh-Ni(100) surface 
using DMol? code (Delley, 1990; 2000). The generalized gradient approximation of Perdew and Wang (GGA-PW91) (Perdew and 
Burke, 1996; Perdew et al., 1996) was used in this study in conjunction with a double numeric quality basis set with polarization 
functions (DNP). The smearing parameter and space cutoff was set to be 0.14 eV and 4.6 A, respectively. Adopting the models of Fan 
et al. (Fan et al., 2012), Li et al. (Li et al., 2014), and Li et al. (Li et al., 2013), asymmetric four-layer slabs, separated by a vacuum of 
10 A, were used in all calculations. However, the binding energies of CH, CH, and H were tested using 4-layer slab fixing the bottom 
two layers and four-layer slab fixing the most bottom layer. The results were similar, where the binding energies using the latter model 
were found to be larger at most by 2.08 kcal/mol. Therefore, the total energies of the slabs were minimized by allowing the top two 
layers of the slab and the adsorbates to move. The tolerances on the energy, gradient and displacement convergence were set to 
2.7 x 10-5 eV, 5.4 x 10-4e V/A, and 5 x 10-3 A, respectively. The dipole slab correction was applied to neutralize the dipole moment 
resulting from using asymmetric slabs. A (3 X 3) unit cells, with a 3 x 3 x 1 Monkhorst—Pack (Monkhorst and Pack, 1976) k-point 
sampling set were used to investigate the reactions and the adsorption of different species associated with them. Transition state structure 
searches were conducted using the Linear (LST) and Quadratic (QST) synchronous transit methods. 


2 Results and Discussion 


2.1 The isolated CH, molecule and clean Rh-Ni(100) surface 


The C-H bond length and the H-C-H angle of an isolated CH4 
molecule were calculated to be 1.099A and 109.471°, (A) (B) 
respectively. The calculated C-H bond length was more 
significant than the experimental value of 1.087 A (Hirota, 1979) 
by 1.1 %. The estimated H-C-H angle is less than the 
experimental value of 109.5° (Van Vleck, 1934) by 0.02%. The 
Rh-Ni(100) surface was modeled by replacing only one Ni atom 
of the most top layer in the unit cell with an Rh atom, as shown 
in Figure 1. The Rh atom was prominent compared to Ni atoms 
in the topmost layer by 0.16 A, as illustrated in Figure 1 (B). The 
average distance between the first and second layers of the slab, 


s . : . Fig.1 (A) Top and (B) side views of a (3 X 3) unit cell of an Rh-decorated 
dı2. was found to be 1.74 A [See Fig. 1 (B)], showing an inward Ni(100) surface. Blue and grey balls represent Ni atoms in the first 


relaxation of 2.7% compared to the bulk distance. The pure and second layers, respectively. The red ball represents the Rh atom. 
Ni(100) surface also shows, using DFT calculations, an inward relaxation of 7.7% (Zhu et al., 2007). Therefore, the presence of Rh 
atom inhibits surface relaxation; however, this could enhance the process of carbon diffusion into the slab, resulting in the formation 
and growth of graphene or carbon nano-tubes (Li et al., 2015a, b). 

Six sites are available on the surface for the adsorption of CH, species as well as carbon and hydrogen, as shown in Figure 2. The 
binding energies between the adsorbates and the adsorbent were calculated using the equation: 
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E, = ~ [Ex/stap = slab —NE,, | (5) 
where Esap, Ex, and Ey/siap are the total energies of a clean slab, an isolated x atom/molecule, and the x/slab adsorption system, 
respectively, with n being the number of adsorbed species within each supercell. The most stable configurations for the chemical species 
on the surface are listed in Table 1 and their binding 


. . j 3 Site Notation 
energies. To determine the reaction energetics, as 
bo Š Bridge between Rh and Ni atoms B_Rh 
well as transition states, the most stable adsorption 
" è $ N y : Bridge between two Ni atoms B_Ni 
configurations for various species associated with the 
CH, dissociation reactions on the surface have been Topon Ratom TBh 
used Top on Ni atom next to Rh atom T Nil 
Top on Ni atom far from Rh atom T_Ni2 
P ä ar Holl Holl 
2.2 CH, adsorption and dissociation said isa 
The following reactions can describe methane Fig. 2 Schematic diagram of a top view of the surface showing the different possible 
adsorption and dissociation adsorption sites on it. Blue and grey balls represent Ni atoms in the first and 
CH ++2 CH (R1) second layers, respectively. The red ball represents the Rh atom. 
4(g) E 4(ads) 
CH. +2 CH +H R2 
acads) 3(ads) i (ads) (R2) Table 1: Most stable adsorption sites for the chemical species on the 
It was found that methane is weakly adsorbed on the Rh-Ni (100) surface and their binding energies in kcal/mol 
surface, where the most stable site is B_Rh (see Fig. Species Most Stable site Binding Energy 
2) with a binding energy of -7.39 kcal/mol; C Hollow -186.19 
inconsistent with the results of Li et al.(Li et al., H Hollow -65.84 
2013) (-6.93 kcal/mol) and Zhu et al.(Zhu et al., cH Hollow “IGLO 
2007) (-2.08 kcal/mol) for CH, adsorpti oe pede apd 
-2. a mol) for CH4 a sorpi ion on CH, B Ni 478) 
NiRh(111) and Ni(100), respectively. Reaction (R2) CH, B_Rh -7.39 


was found to be endothermic by -.32 kcal/mol and 
with a transition state barrier of 21.60 kcal/mol. The reactants, transition state, and products configurations are shown in Figure 3. 
Using NiRh(111), Li et al (Li et al., 2013), and Fan et al. (Fan et al., 2012)found that this reaction is endothermic with a barrier of 17.32 
kcal/mol and 17.79 kcal/mol, respectively. These results indicate that the NiRh(111) surface is more favorable for CH, dissociation. 


2.3 CH} adsorption and dissociation 

The most stable sites for CH, adsorption was predicted to be B_Ni 
followed by B_Rh (see Fig 2) with binding energies of -47.82 and 
46.66 kcal/mol, respectively. The binding energies were 


calculated to be -45.74 and -42.04 kcal/mol for CH, on the Meactants rss neoni 
Ni(100) and NiRh(111) surfaces, respectively. 


The dissociation of CH; is described by the reaction Fig. 3 Initial, transition state and final configurations associated to reaction 
R2. Blue large balls represent Ni atoms, the large red ball represents 


CH3(aasy +*> CHaçaas) + Haas) (R3) Rh atom, the smaller grey ball represents C atom, and white balls 
The reaction was found to be endothermic by 0.76 kcal/mol and represent H atoms. 

has a barrier of 15.5 kcal/mol. The reactants, transition state, and 
products configurations are shown in Figure 4. The reaction was 
reported to be exothermic on the NiRh(1 11) surface with a critical 
barrier of 9.93 kcal/mol (Li et al., 2013). Fan. (Fan et al., 2012) 
found that reaction slightly endothermic with a 13.17 kcal /mol 
barrier energy. This gives the advantage for the NiRh(111) 
surface for this reaction. 


2.4 CH, adsorption and dissociation 


Reactants g Products 
Fig. 4 Initial, transition state and final configurations associated to reaction 
The hollow site was found to be the most stable for CH, R3. 


adsorption with a binding energy of -100.49 kcal/mol. The 
dissociation reaction of CH, is described by the equation: 
CH 2 (ads) +> CH aas) + Haas) (R4) 
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The reaction was exothermic on the Rh-Ni(100) surface by 14.1 kcal/mol and has a 5.96 kcal/mol barrier. The reactants, transition state, 
and products configurations are shown in Figure 5. Fan (Fan et al., 2012) reported that this reaction exhibits an exothermic nature with 
an activation energy of 5.78 kcal/mol. Li et al. (Li et al., 2013) have reported this reaction to be exothermic by 5.47 kcal/mol on 
NiRh(111), however, with a lower barrier of 3.28 kcal/mol, which gives the advantage to the (111) configuration for CH, dissociation. 


2.5 CH adsorption and dissociation 


The hollow site was the most stable site for CH adsorption, with 
a binding energy of -161.70 kcal/mol. The CH dissociation takes 
place according to the reaction 

CH (aas) FAD C(aas) + Heaas) (R5) 


The calculations show that this reaction has an exothermicity of 


Reactants TS3 Products 


1.22 kcal/mol and requires a barrier of 16.8 kcal/mol. The Fig. 5 Initial, transition state and final configurations associated to reaction 
reactants, transition state, and products configurations are shown R4. 

in Figure 6. Reaction (R5) was stated to be endothermic on 

NiRh(111) by 15.25 kcal/mol and requires a barrier of 29.8 kcal/mol (Li et al., 2013). Fan (Fan et al., 2012) have also found that this 
reaction is endothermic and requires a 32.11 kcal/mol activation energy. These results show that the Rh-Ni(100) surface is more adroit 
for the CH dissociation reaction. Atomic carbon and hydrogen were found to be adsorbed on the hollow site with binding energies of 


-186.19 and -65.84 kcal/mol, respectively. 


Reactants TS4 Products 


2.6 The overall reaction 


Figure 7 shows the energies along with the coordinates of the 
overall reaction CHy(g) = Caas) + 4H(aas) over the Rh-Ni(100) 
and NiRH(111) surfaces. Figure 7 gives a better understanding 
of the whole process and explains which surface best facilitates 


the whole CH, dissociation process. : n N i : ; : ; 

i j 2 Fig. 6 Initial, transition state and final configurations associated to reaction 
From the discussion above, the NiRh(111) surface was found to R5. 
be more favorable for the reactions (R2), (R3), and (R4), whereas 
Rh-Ni(100) was found more preferable for reaction (R5) for the lower required barriers. This is true if these reactions were carried out 
individually, but the reactions simultaneously occur on the surface. The complete process, therefore, should be taken into account, and 


here, two particular processes will be discussed; these are CH} and CH dissociation reactions. The dissociation of CH3 can be described 


Eq_ Eq 
by the equation CH, 2 CH, + H ms CH, + 2H, where Eap; is the activation energy of the forward reaction of R3 and E,_,, is the 
activation energy of the backward reaction of R2. On the NiRh(111) surface, E,,, was found to be 9.93 kcal/mol (Li et al., 2013), lower 
than the activation energy required for the same reaction on the Rh-Ni(100) surface by 5.57 kcal/mol. However, E,_,, for the 


reformation of CH, was found Rh-Ni(100) NiRh(111) 
. 30 30 (23.79) 
to be 9.70 kcal/mol using S4 
‘ i . 20! 20l (16.63) (16.86) 
NiRh(111), which is lower p phere (12.43) _ (st) (182) 
than the barrier needed for £ 10 i — £ 10 
: or = i = 9.24) 
CH; dissociation by 0.23 § d pel o is ey = oan) | 
kcal/mol. This means that if > ` (9.60) (CH,+2H) 
š D D 4 
a M (-6.01) 
there is a hydrogen atom next 2 10 E 10 Pr 
to the CH, entity on the (4641) 
: -20r Cham)  €1763)] -20 
NiRh(111) surface, the (C+4H) 
reformation of CH, will be oe Reaction Coordinates = Reaction Coordinates 
more likely to occur. On the Fig. 7 The potential energy surfaces for CH, dissociation on Rh-Ni(100) and NiRh(111) surfaces. 


other hand, the barrier of CH, 
reformation reaction on the 
Rh-Ni(100) surface is calculated to be 17.28 kcal/mol, which is higher than the barrier of CH3 dissociation by 1.78 kcal/mol, giving the 
advantage to the reaction to go further towards CH, and H production. Thus, the Rh-Ni(100) surface is anticipated to be better for this 
reaction than NiRh(111), notably, if an atomic hydrogen coexists with CH3 on the surface. The dissociation of CH was reported to be a 
limiting step for CH, decomposition reaction on NiRh(111) due to the required activation energy (Li et al., 2013; Zhu et al., 2007). The 
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a-R4 


E Ea 
process of CH decomposition can be described by: CH, = CH + H 3 C + 2H where E 
reaction of R5 and E,_,, is the activation energy of the backward reaction of R4. 

The value of Eag, was calculated to be 16.80 kcal/mol and 29.78 kcal/mol using the Rh-Ni(100) and NiRh(111) surface, respectively, 


giving the advantage for the Rh-Ni(100) surface. Besides, the dissociation of CH on NiRh(111) requires a barrier of 29.8 kcal/mol, 
whereas the barrier of the backward reaction of CHzçaas) +*= CH (aas) + Hcaas) needs a barrier of 12.70 kcal/mol (Li et al., 2013), 


aps 18 the activation energy of the forward 


which means that if atomic hydrogen coexists next to an adsorbed CH entity, the reformation of CH, is more likely to take place instead 
of CH splitting. Using Rh-Ni(100), the barrier of CH dissociation is lower than the activation energy required for CH, reformation by 
16.3%, giving the benefit to the reaction to go further towards C and H formation, leading therefore, to a higher conversion of CH to C 
and H. 


Conclusions 


The mechanism of the reaction CH4(g) = C(aas) + 4H(aas) Was investigated using DFT calculations for the periodic slab on the Rh- 
Ni(100) surface. It was shown that methane is weakly adsorbed on the surface and requires a 21.60 kcal/mol barrier to dissociate. The 
results of this investigation were compared to the results of the same reaction over NiRh(111). It was shown that the NiRh(111) 
configuration is more favorable for reactions R2, R3, and R4 if each reaction was carried out individually on the surface, whereas Rh- 
Ni(100) is better for the reaction R5. However, the Rh-Ni(100) surface was found to benefit the complete process in a better way, due 
to the decrease in the energy level of barriers along with the reaction coordinates, which leads to a higher conversion of CH, to carbon 
and hydrogen. 
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Nomenclature 

Egr =the activation energy of the backward reaction [kcal/mol] 

E, =the activation energy of the forward reaction [kcal/mol] 

Estab =total energy of a clean slab [kcal/mol] 

Ey =total energy of isolated x atom/molecule [kcal/mol] 

Ex/siap =total energy of the x/slab adsorption system [kcal/mol] 

n =the number of adsorbed species within each supercell [atom or molecule] 
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